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ABSTRACT 

We report on a 10 ksec observation of the double pulsar system J0737-3039 obtained with the Chandra 
X-ray Observatory's Advanced CCD Imaging Spectrometer. We detect the system as a point source with 
coordinates (J2000) RA = 07'^37'"5L;23, Decl. = -30°39'40'.'9, making this the first double neutron star 
V-j.". system to be detected at X-ray energies. Only 77 source counts are detected. The data are represented 

' satisfactorily by a power-law spectrum with photon index F = 2.9 ± 0.4. The inferred 0.2-10 keV 

I luminosity is « 2 x lO'^^erg s~'^ for a distance of 0.5 kpc. We do not detect any significant variability 

CS| ■ with orbital phase in these data, which span one full orbit. These results are consistent with the X- 

ray emission originating solely from the magnetosphere of the energetic A pulsar, with an efficiency for 
, conversion of rotational energy to X-rays in this band of ^ 0.04%, although we cannot exclude other 

[t I ■ possibilities. 

^ ■ Subject headings: pulsars: individual (PSR J0737-3039A, PSR J0737-3039B) — stars: neutron — 

^SJ I X-rays: stars 



1. INTRODUCTION 

Thirty years following the discovery of the first double 
neutron star (Hulse & Taylor 1975), only eight such sys- 
tems are known. A few of these rare binaries are wonderful 
laboratories for the study of relativistic gravity and gravi- 
tational radiation. One of the more recently discovered of 
these systems, J0737-3039 (Burgay et al. 2003), is funda- 
mentally different from all of the others observationally in 
that it contains two detected pulsars (Lyne et al. 2004). 
The first pulsar discovered in the system (A) spins every 
22 ms, has a low inferred surface magnetic field strength, 
a characteristic age of 210 Myr, and a relatively high rate 
of rotational energy loss, E = 5.8 x lO'^'^erg s^^. Pulsar 
B has a period of 2.7 s, a characteristic age of 50 Myr, and 
an E that is 3600 times smaller than that of A. Its mass, 
1.25 Mq, is the lowest known for any neutron star and is 
0.09 Mq less than that of A. Both pulsars orbit their com- 
mon center of mass every 2.4 hr with eccentricity of 0.09. 
We observe the system nearly edge-on, with an orbital in- 
clination of 87°. 

A unique aspect of the system is that radio emission 
from pulsar B is only detected strongly in two repeatable 
orbital phase ranges, each of duration lOmin. A 30 s- 
long eclipse of pulsar A is also observed at its superior con- 
junction. Both of these effects most likely result from the 
interaction of the fast-spinning pulsar's relativistic wind 
with the magnetosphere of its much less energetic com- 
panion (Lyne et al. 2004). 

This system will be an exceptional laboratory for 
relativistic astrophysics and should provide unique in- 
sights into the magnetospheres and the local environ- 



ment of both neutron stars. High-energy observations, 
particularly at X-ray energies, could be especially valu- 
able for understanding the energetics of the system. 
Whether resulting from magnetospheric or thermal pro- 
cesses, most pulsars emit of order 10^^ E in soft X-rays 
(Becker & Triimper 1999). This translates to an expected 
luminosity of order lO'^^ erg s~^ for pulsar A, and means 
that B should be undetectable due to its very low E and 
relatively large age. But in this particular system, while 
still unlikely, it is conceivable that a portion of A's rela- 
tivistic Poynting/particle flux might heat the surface of B, 
perhaps causing it to emit detectable radiation. Such emis- 
sion could prove invaluable for probing the atmosphere of 
B, which is likely to deviate significantly from standard 
models. A more likely source of "unusual" X-ray emission 
might be the shock that undoubtedly forms near pulsar 
B owing to the collision between A's relativistic wind and 
B's magnetosphere. Some of this emission could be time 
variable, for instance due to Doppler boosting. 

Many of these emission modes are clearly speculative 
and hard to quantify with the present understanding of 
the system. Nevertheless, given the potential for un- 
usual insights into the physics of this unique interacting 
binary, we have obtained a short X-ray observation of 
this system, awarded through the Chandra Director's Dis- 
cretionary Time program, in order to better plan future 
multi-wavelength studies. In this Letter we report the de- 
tection of a point source at the position of the J0737-3039 
system and describe its properties. 
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2. X-RAY OBSERVATIONS AND ANALYSIS 

The double pulsar system J0737-3039 was observed with 
the Advanced CCD Imaging Spectrometer S-array (ACIS- 
S) instrument aboard Chandra on 2004 January 18 for 
lOksec (corrected exposure time was 10,009 s). The sys- 
tem's position was located on the aim point of the ACIS 
S3 chip. We used the CIAO package^, version 3.0.2, to per- 
form much of the analysis of these data, using the most re- 
cent release of the calibration software, CALDB 2 . 26, which 
corrects for the loss of quantum efficiency at low ener- 
gies due to molecular contamination. As shown in Fig- 
ure 1, we detect a source at (J2000) RA = 07'^37'"51!23, 
Decl. = -30°39'40'.'9. No other sources are detected in the 
entire S3 chip. This position is C.'S away from that derived 
from radio timing of the pulsars (Lyne et al. 2004), well 
within the X-ray 1 a position imcertainty of O'.'S, domi- 
nated by the satellite pointing error of approximately 0"7. 
We have therefore detected the J0737-3039 system in X- 
rays. Using a circular aperture of radius 2" centered on 
the position of the system, and four background regions, 
we determine that 77 ± 9 counts originate from the pul- 
sar, corresponding to a rate of (7.7 ±0.9) x 10""^ counts/s. 
The system is detected as a point source, with the spatial 
spread consistent with the point spread function of the 
instrument. 

We used psextract to obtain a spectrum for this source 
and fit it using Sherpa, the CIAO modeling and fitting pack- 
age. As shown in Figure 2, a steep power-law spectrum 
fits the data satisfactorily, with photon index F 2.9 and 
neutral hydrogen column density Nh « 5 x lO^^cm"^. 
Because of the small number of counts we use the maxium 
likelihood method to fit these data, implemented in Sherpa 
with the CSTAT statistic (Cash 1979). This fit returns a 
Q- value (i.e., the probability that we would observe the 
measured CSTAT statistic or higher if the best-fit model is 
true) of 0.96. The best-fit blackbody and Raymond-Smith 
models return Q- values of 3 x 10~^ and 2 x 10~^. A thermal 
bremsstrahlung model returns a Q-value of 0.55, but the 
best-fit value of Nh ^ 1 x 10^^ cm^^ is unreasonably low. 
Fitting for multi-component models or models with more 
free parameters is unreasonable due to the small number 
of counts available. While there are only two counts with 
energies greater than 3 keV, the fitted parameters do not 
change substantially if these counts are removed. Simi- 
larly, the fit does not change substantially if we exclude 
photons with energies less than 0.5 keV, where ACIS-S 
calibration is somewhat uncertain. 

Figure 3 shows the a contour plot of the CSTAT statistic 
for the best-fit model in the Nh — F plane. From the de- 
rived count rate and spectral fit, we calculate a 0.2-10 keV 
unabsorbed flux of w 8 x 10~^^erg cm~^ s~^. As- 
suming isotropic emission, and for a distance d ^ 0.5 kpc, 
inferred from its dispersion measure (i.e. free electron col- 
umn density) and a model for the Galactic electron density 
(Cordes & Lazio 2002), we obtain a 0.2-10 keV luminosity 
of « 2 X 103°((i/0.5 kpc)2 erg s'^ In Table 1 we list all 
relevant spectral quantities and their uncertainties. We 
stress that d, and hence the derived luminosity, may be 
uncertain by a factor of two or more due to unmodeled 
electron density features associated with its position near 

* Available at http://asc.harvard.edu/ciao. 

^ See the ISM Column Density Search Tool at http://archive.stsci.ei 



the Gum Nebula (Mitra & Ramachandran 2001). VLBI 
observations will be important for determining the actual 
distance. 

The value of Nh derived for the system, while rather 
uncertain, is much smaller than the total Galactic value 
of 4.65 X 10^^ cm~^ derived by Dickey & Lockman (1990) 
from neutral hydrogen measurements. It is also smaller (at 
> la) than the value of 1.5 x 10^^ cm~^ derived from its 
dispersion measure, assuming 10 H atoms for each e~ . Our 
derived Nh. however, is marginally greater' than expected 
given the measured neutral hydrogen absorption to other 
sources with similar longitude, latitude and ~ 0.5 kpc dis- 
tance®, and is more consistent with the measurements if 
a distance of at least ^ 1 kpc is assumed. Such discrep- 
ancies between Nh measurements and dispersion-measure 
derived distances for pulsars are not unusual (e.g. Nicastro 
et al. 2004). Also, since the source closest to J0737-3039 
with a neutral hydrogen absorption measurement is a ra- 
dial distance of ~ 0.1 kpc away, and because the errors on 
our Nh and distance measurements are large, further ob- 
servations are necessary to determine if such a discrepancy 
exists. 

In Figure 4 we plot the number of photons detected 
in each of 10 time bins over the lOksec observation (i.e. 
1.13 orbital periods). There is no significant evidence for 
variability, with a x test revealing only a 51% probability 
that the data were drawn from a non-uniform distribution. 
We also do not sec any correlation between time and the 
energy of the photons. The small number of detected pho- 
tons makes searching for variability on shorter time scales 
impossible. 

3. DISCUSSION 

We have detected the J0737 3039 system as a faint X- 
ray source that is best described by a steep power-law spec- 
trum with F « 2.9 and » 2 x 1030(d/0.5kpc)2erg s'^ 
(see Table 1). This result is compatible with the emis- 
sion originating solely from the magnetosphere of pulsar 
A: Becker & Triimper (1999) find that for recycled radio 
pulsars that are detected at X -ray cncrgiGS, J-^x 
and that when their spectra can be fit by a power law, 
2 < F < 2.4. While by comparison our spectral index is a 
bit steep and our efficiency for conversion of rotational to 
X-ray energies is a little small (4 x 10~"^[d/0.5kpc]^), un- 
certainties in spectral parameters and distance make this 
interpretation perfectly plausible. We note in passing that 
the X-ray luminosity detected from the system is roughly 
equal to the entire spin-down luminosity of pulsar B. 

Another intriguing possibility is that these X-rays are 
produced when the relativistic winds of A and B collide. 
As shown by Arons & Tavani (1993), the synchrotron emis- 
sion produced in such interactions would also be expected 
to have a power law X-ray spectrum with F ~ 2. Un- 
fortunately, the data from this short observation are not 
sufficient to discriminate between these two possibilities. 
More sensitive observations with better time resolution are 
necessary in order to: obtain a better spectral fit; deter- 
mine what portion of the detected X-ray flux is pulsed at 
the rotation period of pulsar A; quantify better any pos- 
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sible time variability (the observed radio flux density of 
B varies dramatically on ~ 5-10 min timescales that are 
impossible to investigate with the present data); measure 
any thermal component that may be contributing to the 
total X-ray flux of this system. 
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Table 1 

x-ray properties of the j0737-3039 system. 



Parameter Value 

r, photon index '^-^to.t 

Nfj, neutral hydrogen column density (10^°cm~^) 4.812 4 

fa:, flux (0.2-lOkeV) (IQ-i'^erg cm-^g-i) 7.6tii 

La:, luminosity (0.2-10 keV) (10^° erg s"^) 2.3t^ j 

La:/E (0.2-10 keV) (lO"-*) 4.0t^-^ 



Note. — Results from a fit to ACIS-S3 data of the form 
dN DC E~^dE. Flux has been corrected for interstellar absorp- 
tion. Luminosity assumes a distance of 0.5 kpc (see § 2) and is 
the same, within the errors, for the 0.2 - 3 keV energy range. 
Uncertainties correspond to 1 ct confidence levels. 
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Fig. 1. — A 50" X 50" portion of the ACIS-S3 chip in the 0.2-10 keV band, with the double pulsar J0737-3039 at its center. No other 
sources were detected in the entire 16' X 16' area of the chip. 
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Fig. 2. — The background-subtracted spectrum of J0737-3039, with the power-law fit described in § 2. Each of the first five energy bins 
contains 15 counts; the last bin contains two counts. The squares are placed at the average energy of the photons in each bin, while the 
solid best-fit line illustrates each bin's full width. Error bars, calculated by Sherpa, account for uncertainties in both source counts and the 
instrumental model. 
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Fig. 3. — Contour plot of the CSTAT statistic as a function of the neutral hydrogen column density Nh and the photon index F. The 
plotted confidence levels are 1, 2 and 3it. 



-i-J 

O 
O 




2000 



4000 



Time (seconds^ 



Fig. 4. — Number of counts (shown with Poisson error bars) vs. elapsed time (since the start of the observation) with 10 bins across 
the 10 ksee observation. The orbital period of the system is 8835 s. Orbital phase calculated from barycentercd data (with respect to the 
ascending node of pulsar A) ranged from 281° at the start of the observation to 328° at the end, with the dashed line marking the superior 
conjunction of A. 



